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Abstract 

High Mass X-ray Binary Pulsars (HMXBP), in which the companion star is a 
source of supersonic stellar wind, provide a laboratory to probe the velocity and 
density profile of such winds. Here, we have measured the variation of the absorption 
column density along with other spectral parameters over the binary orbit for two 
HMXBP in elliptical orbits, as observed with the Rossi X-ray Timing Explorer 
(RXTE) and the BeppoSAX satellites. A spherically symmetric wind profile was 
used as a model to compare the observed column density variations. In 4U 1538-52, 
we find the model corroborating the observations; whereas in GX 301-2, the stellar 
wind appears to be very clumpy and a smooth symmetric wind model seems to be 
inadequate in explaining the variation in column density. 
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4U 1538-52 (Giacconi et al. 1974) and GX 301-2 (White et al. 1976) are both 
luminous (~10^®-10^^ ergs s“^) pulsars in elliptical orbits. They emit X-rays 
due to the accreted matter from the stellar wind of their respective companion 
stars (Parkes et al. 1978, Kaper et al. 1995). Since both have well dehned 
orbital measurements (Glark 2000, Koh et al. 1997), it enables us to inspect 
the variation of the X-ray spectral parameters at different orbital phases to 
gain useful insight into the morphology of the stellar wind. In this work, we 
have measured the X-ray spectral evolution of these two HMXBP as observed 
with the RXTE and BeppoSAX satellites and the parameters of the spectral 
model were examined along the orbit for any possible variations. We have also 
compared the equivalent hydrogen column density (Nh) variations over the 
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orbit with a model variation, assuming a spherically symmetric stellar wind 
emanating from the companion star. 

We observed 4U 1538-52 with RXTE from 2003-07-31 to 2003-08-07 covering 
out of eclipse phases for two binary orbits. We also used the archival data from 
BeppoSAX obtained between 1998-07-29 to 1998-08-01, covering one binary 
orbit. GX 301-2 was observed by RXTE, Erst from 1996-05-10 to 1996-06-15 
and again from 2000-10-12 to 2000-11-19. For more details, refer to Mukherjee 
and Paul (2004) and Mukherjee et al. (2005) 

For the RXTE data, we took Standard 2 data products of the Proportional 
Counter Array (PGA) and extracted the source spectra using the tool saex- 
trct V 4.2d. The BeppoSAX data products were extracted from the Medium 
Energy Concentrator Spectrometers (MEGS) and Low Energy Concentrator 
Spectrometers (LEGS) using circular regions of radii 4' and 8' respectively. 
The background subtracted source spectra were analyzed with XSPEC v 
11.2.0 (Shafer, Haberl & Arnaud 1989) using appropriate spectral models. 
For 4U 1538-52, we used a power law along with a line-of-sight absorption and 
a gaussian line with centre energy ~ 6.4 keV. We also applied a high energy 
cut-off component. And for GX 301-2, we used the Partial Covering Absorber 
Model (PCAM, Endo et al. 2002) along with the addition of a high energy 
exponential cutoff. The PCAM is described as two different power law com¬ 
ponents with the same photon index but different normalizations (Norml & 
Norm2); being absorbed by different column densities, Nhi & Nh 2 respectively. 
Xh 2 is interpreted as the column density of the material local to the X-ray 
source, while Nm accounts for the rest of the material (along with galactic 
absorption). In this case, iron Ka and K/J lines and an absorption edge at 7.1 
keV due to neutral iron were also included in the model. 


Results &; Discussion 


4U 1538-52 : The photon index, fluorescent iron-line flux, cut-off energy and 
the e-folding energy measured with the average spectrum taken over 2-3 ks 
does not show any substantial variation along the orbit. This suggests that the 
continuum X-ray spectrum of the pulsar is hardly affected during its revolution 
along the orbit. But we detected a notably smooth variation in Nh with orbital 
phase, increasing gradually by an order of magnitude as the pulsar approaches 
eclipse (near phases 0.3 and 0.5; as measured from the periastron passage. Fig. 
1). At orbital phases far from the eclipse, the column density has a value of 
~ 1.5 X 10^^ H atoms cm“^. We compare the observed column density profile 
with a model estimated by assuming a spherically symmetric Gastor, Abbott 
& Klein (GAK 1975) wind from the companion star. The velocity profile of the 
wind is : v^ind = where Voo : terminal velocity of the stellar wind. 
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Fig. 1. The variation of column density (10^^ H atoms cm“^) with orbital (binary) 
phase for 4U 1538-52 (left) and GX 301-2 (right) are shown. In the former, the 
dashed line, the dashed-dotted and the dotted line represents the model Nh values 
for the inclination angles of 65°, 75° &: 85° respectively. In the latter, the three 
panels show the model variation due to the wind and the measured variations of 
Nhi and Nh 2 respectively. In both the figures, different symbols are used for different 
observations. 


Rc : radius of the companion and r : radial distance from the centre of the 
companion star. The column density prohles were derived using a numerical 
integration along the line of sight from the pulsar to the observer for three 
different inclination angles of 65°, 75° and 85° respectively. With a mass-loss 
rate of Mq yr“^ and Vqo ~1000 km s“^, the model calculations of Nh for 

different inclination angles when superposed on the observed values show fairly 
reasonable agreement (Fig. 1), indicating that a spherically symmetric stellar 
wind from the companion star may produce the observed orbital dependence 
of the column density for certain range of the orbital inclination. 


GX 301-2 : In most cases, the spectrum was htted well with the PCAM. 
The average values of the free parameters measured here over the full binary 
orbit; viz. photon index, e-folding energy and cut-off energy follow the general 
trend which were earlier measured only in some phases of the binary period 
(White et al. 1983, Orlandini et ah 2000). But the variation of Nhi & Nh 2 with 
orbital phase was not smooth. The values were very high with a large varia¬ 
tion throughout the binary orbit (from 10^^ to 10^"^ atoms cm“^), indicating a 
clumpy nature of the stellar wind (Fig. 1). It is also seen from Fig. 2 that the 
covering fraction (C.F., dehned as Norm2/[Norml-|-Norm2]) remains substan¬ 
tially high almost throughout the orbit. This means that there is dense and 
clumpy material present all through. As in the case of 4U 1538-52, we com¬ 
pared the measured values of column densities with a model variation using 
a CAK wind. The model values were of the order of 10^^ to 10^^ atoms cm“^ 
(Fig. 1). The peak between phases 0.1 and 0.2 was expected since the line of 
sight passes through the densest parts of the wind during these phases. Thus 
it is clear that the observed variation in column density cannot be explained 
by a spherically symmetric CAK wind only, indicating probably strong inho- 
mogeneties in the wind. Now, as can be seen from Fig. 2, the two iron lines 
included in our model show large increase in flux near periastron (phase ~ 0.9) 
and a possible small increase near phase 0.1 (at least for 1996, hlled squares). 


3 










0.03 




4'_A -a ■ # - a ■ ■ ■ -A A^^ 

A ^ A ■ 

0.02 

6.4 keV Flux 

A 

4* ' 

A A ^ - 

& . ^ 

n" 0.01 


4 

4- 

A ^ 

o 0 

_ a' tR* A- ^ . 4. ilf . i ^ 

. 4 ^ 

& 

A A 

6 xl 0 “® 


4 

A 

i 

ft; 

A ■ 

4x10 ^ 

3x10“® 

7.0 keV Flux 

-i* J A* 14 ** 

.4, t 

0 0.2 0.4 0.6 0.8 1 

0 

4 


Binary Phase 0 0.3 0.4 0.6 0.9 1 

Binary Phase 


Fig. 2. The variation of C.F. (left) and the iron-line flux (right) with the binary 
phase for GX 301-2 are shown with different symbols for different observations. 

The peak near periastron is not very evident in the 1996 data due to the lack 
of enough observations, though an increasing trend can possibly be inferred. 
Moreover, the equivalent width of the 6.4 keV iron-line showed a correlation 
with Nm, suggesting that most of it is produced by the local clumpy matter 
surrounding the neutron star. 
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